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We investigate decoherence of an electron in graphene caused by electron-flexural phonon inter¬ 
action. We find out that flexural phonons can produce dephasing rate comparable to the electron- 
electron one. The problem appears to be quite special because there is a large interval of temperature 
where the dephasing induced by phonons can not be obtain using the golden rule. We evaluate this 
rate for a wide range of density (n) and temperature (T) and determine several asymptotic regions 
with temperature dependence crossing over from ~ to ~ T when temperature increases. 
We also find to be a non-monotonous function of n. These distinctive features of the new contri¬ 
bution can provide an effective way to identify flexural phonons in graphene through the electronic 
transport by measuring the weak localization corrections in magnetoresistance. 

PACS numbers: 72.10.-d, 72.10.Di, 72.80.Vp 


Introduction. The transport properties of graphene 
have attracted much attention [1] since the first discov¬ 
ery of this fascinating material [2]. It is promising for 
various applications due to its high charge mobility and 
unique heat conductivity. Theoretically, it was realized 
long ago [3-5] that these transport properties of free¬ 
standing (suspended) graphene are strongly influenced 
by flexural (out-of-plane) vibrational modes that deform 
the graphene sheet. From the experimental point of view, 
the effect of flexural phonons (FPs) was clearly observed 
in heat transport [6, 7]. However, it is a more challeng¬ 
ing task to identify the effect of flexural phonons in elec¬ 
tronic transport [8, 9]. This is because the contribution 
of electron-phonon interactions to momentum relaxation 
remains small even at high temperatures, with the main 
source of the relaxation being elastic impurities [10]. 

The dephasing rate on the other hand, is a more 
suitable quantity for studying FPs, since static impuri¬ 
ties do not cause dephasing. Usually, electron-electron 
interactions, [11-15] are considered the primary mecha¬ 
nism for dephasing. In this letter we discuss dephasing 
caused by the electron-flexural phonon (el-FP) interac¬ 
tion in graphene. It is the softness of the flexural mode 
and the coupling of an electron to two FPs simultane¬ 
ously (see Fig. 1 for illustration) that make the contri¬ 
bution of FPs to significant in a suspended sample, 
and at large enough densities comparable with the one 
caused by the electron-electron interaction. Because of 
the quadratic spectrum of FPs, ujk = ak^, they are much 
more populated as compared with in-plane phonons. In 
addition, the coupling to two FPs considerably increases 
the phase space available for inelastic processes as com¬ 
pared to the interaction with a single phonon. The point 
is that in graphene the Fermi momentum, kF,is relatively 
small. As a result, the interaction of a single phonon with 



Figure 1. On the left: scheme of the el-FP interaction process, 
where the solid line represents an electron, and the wavy lines 
represent FPs. On the right: FPs can have momenta p, q 
much larger than the transferred momentum Q. Under the 
conditions discussed in the paper, the scattering process is 
considered as semi-elastic. 


electrons is determined by the Bloch-Griineisen temper¬ 
ature, Tbg ^ 2 kpi rather than the temperature, when 
T ^ Tbg [16]- In such a case, one needs to exploit other 
scattering mechanisms to overcome the limitations in¬ 
duced by the smallness of kp [17]. In the case of el-FP in¬ 
teraction, coupling to two phonons radically changes the 
situation. Now only the transferred momentum should 
be small, while individually a FP may have a momentum 
much larger than kp, up to the thermal momentum qp. 

Still, as we shall demonstrate, the problem of dephas¬ 
ing due to the el-FP interaction appears to be quite spe¬ 
cial, because the softness of FPs, i.e. unique smallness of 
TbG) leads to the existence of a temperature range where 
dephasing rate cannot be obtained using the golden rule 
(GR). Rather, both the self-energy and the vertex pro¬ 
cesses [18] should be treated simultaneously. This results 
in a transition from ~ to T with increasing tem¬ 
perature for the dephasing rate induced by FPs. 

The electron-flexural phonon interaction. Lattice dy¬ 
namics of the single-layer graphene can be described 
in terms of the displacement vector u = {ux,Uy,h) [19]. 




2 


Here describe the in-plane modes, while the out-of- 
plane displacement h describes the flexural mode. The 
displacement vector leads to a non-linear strain tensor 
Uij = ^ {diUj -I- djUi -I- dihdjh ), where (*,j) = {x^y) are 
spatial indices. The lattice modes interact with elec¬ 
trons through emergent scalar and vector potential fields 
[ 20 , 21 ]: 


^ - Ql T Uyy') , 

A — S g2 / Xp {Uxx ^yy ? ) , (1) 

where gi = 30eV, 52 = 7.5eV [9] and is the Fermi 
velocity. Index a = K,K' describes two valleys of the 
conducting electron band, and factor = ±1 reflects 

the fact that the emergent vector potential A“ respects 
the time reversal symmetry. 

Thermal fluctuations of the lattice produce variations 
in the potentials. Averaging over lattice vibrations one 
finds the correlation functions of the potentials as 

((p(Q,H)(/?(-Q,-H)) =(/)(Q,H), 

(a“ (Q, H) Af (-Q, -H)) = (Q, H). (2) 

To proceed, we introduce the correlation function for FP 

{h (k, w) h (— k, —uj)) = H ( fc ) 2tt5 {uj — Uk) , (3) 

where H [k) = In this equation, n (w) is the 

Planck distribution function and p is the mass density 
of the graphene sheet. One can propose the following 
form of the spectrum of the flexural phonon: 

ccfc = ak^Q {k), Qik) = ^l + Z-^{qJk)\ (4) 

where 0 (fc) describes a transition from the bare spec¬ 
trum at high momentum to the renormalized spectrum 
~ /;;2-7)/2 momentum limit. At k<qc (T) = 

the quadratic spectrum for the flexural mode 
ceases to work due to anharmonicity. Here Ac « 18.7eV 
[5] reflects the energy scale of anharmonicity. The an¬ 
harmonicity is related to the /I'^-vertex, arising as a re¬ 
sult of integrating out fast u-modes, which are coupled 
to /i-mode [22]. Below we will exploit the value Z ~ 2, 
and take 77 « 0.8 from the numerical solution of the self- 
consistent screening approximation theory [23, 24]. 

We consider graphene away from the Dirac point at 
chemical potential p, ^ T. Besides kp, the relevant 
momentum scales in the problem are: thermal mo¬ 
mentum qx = i/r/a « 0.05\/T[K]/A, and qdT) « 
0.01-\/T[K]/ a which signals the transition to the renor¬ 
malized FP spectrum. From now on, we will concentrate 
on the realistic situation from the experimental view¬ 
point: kp ^ qx, i.e., T 3> Tpa- The Bloch-Griineisen 
temperature Tbg = ^ 2 kF ~ 0.40 {2kp) nK, where n is 
the electronic density measured in units of 10^^cm“^. 


Note Tbq is extraordinarily small for all relevant den¬ 
sities. As we have already emphasized, see Fig. 1, the 
momentum transfer in the el-FP interaction is limited 
by 2kp. Nevertheless, the extended structure of the cor¬ 
relation functions (j) (Q, H) and Aij (Q, H) enables elec¬ 
trons to have energy transfer exceeding the phonon en¬ 
ergy UJ 2 kF ■ 

Main tool to probe electronic coherence is magnetore¬ 
sistance [25], which gives a direct access to the weak lo¬ 
calization corrections to conductivity, controlled by the 
dephasing rate . The weak localization correction to 
conductivity in graphene can be written as [26, 27] 

Aa = - — ^ / dt& {-t/2, t/2 ), (5) 

^ i '' 

where I sums over four Cooperon channels relevant for the 
magnetoresistance. Physically, C* (—t/2, t/2) represents 
the interference of a pair of time reversed trajectories in 
the channel I that start at —t/2 and return to the ini¬ 
tial point at t/2. More generally, the Cooperon matrix 

labelled by two isospin numbers si ,2 and two 
pseudospin numbers li^ 2 - This matrix is diagonal in the 
pseudospin space even in the presence of interactions that 
preserve sublattice and valley indices. The Cooperon 
channels relevant for magnetoresistance are the isospin 
singlets, C* = Cqq (Z = 0,x,j/, 2 :), that do not have gaps 
comparable with r“^, the elastic scattering rate due to 
impurities. Therefore, we restrict ourselves to this sub¬ 
space. 


(a) 


(b) 



Figure 2. Diagrammatic representation of (a) the self-energy 
and (b) the vertex FP-contribution to the Cooperon; see also 
Fig. 2c in [28]. 

To include el-FP interaction into the Cooperon, one 
can write down a Bethe-Salpeter equation for a particular 
Cooperon channel C^, see Fig. 2. In the following we 
will not solve the equation exactly, but instead, we will 
estimate the upper bound of the Cooperon decay rate 
[29, 30]. We start by writing down an ansatz that reads 
as [18] 


( 6 ) 
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Here Cg (t) is the diffusion propagator describing the bare 
Cooperon, and F* (^ 1 ,^ 2 ) is a decay function characteriz¬ 
ing the effect of the el-FP interaction [28]. 

Dephasing due to scalar potential fluctuations. For the 
scalar potential correlation function one obtains: 

(f) (Q, H) =^gl (Q) J (d^p) jp . ^^2 

xF(p)F(g)dp,q(H,Q), (7) 

where dp,q (H, Q) = (27r)^ d (H ± Wp ± Wq) x 

(5(Q —p —q), and Wp^q are given by Eq. (4). Here 
summation includes four different processes of emis¬ 
sion/absorption of two FPs by an electron. The screened 
coupling constant gi {Q) = 51 where x = g^Nkp, 
IV = 4 is the spin-valley degeneracy in graphene, and 
ge 1 describes the renormalized Coulomb interaction 
[31]. Since each time an electron is coupled to two 
flexural phonons, (j) describes a phonon loop and, there¬ 
fore, in the momentum-frequency domain (j) (Q, H) has a 
extended support rather than a d-function peak. As a 
result, the decay function for the scalar potential (t) 
(which is the same for all channels) can be expressed as 
a convolution of the three factors [28]: i) the correlation 
function ^(Q, H), ii) function describing the 

ballistic electron’s motion, and hi) factor {fl,t), 
reflecting the relation between the self-energy and vertex 
diagrams: 

it)=t I (dQ) (dH) <). (Q, n) (Q) (H, t). ( 8 ) 

Here, 

(Q) = ^ (1 - 0 (2fcF - Q) , (9) 

where the Heaviside theta function 6 (2kp — Q) restricts 
momentum that can be exchanged between FPs and elec¬ 
trons. The factor (H, t) is equal to 

C^( 0 ,i) = l-^, ( 10 ) 

and it describes the balance between the self-energy and 
vertex diagrams on Fig. 2. is sensitive to dynamic 
aspect of the scattering event and, because of this, alters 
temperature dependence of . 

The dephasing rate is defined according to 

F<l>(T(f,) = 1- The decay function can be most conve¬ 
niently expressed as 

F^it) = cltTfir,o-, ( 11 ) 

where = ^ 2 -Kg°N 1-2 1 ® dimensionless coupling con¬ 

stant and / is a dimensionless function of two parame¬ 
ters: T = akpt and ^ = Z~^l^qc/kF [28]. Parameter ^ 


T(A-) 



Figure 3. Phase diagram of the dephasing rate due to FPs 
with scalar coupling. The blue and black lines divide the 
whole (T — fj,) plane into three regions, see the text for ex¬ 
planations. The blue line coincides with the maximum of the 
dephasing rate as a function of chemical potential at a fixed 
temperature, see Fig. (4). The red dashed line representing a 
fragment of ^ = 1 is shown here for orientation. The inset is 
a zoom in of the intersection area of the blue and black lines 
plotted as a function of the electronic density. 

originates from the renormalization of the FP spectrum 
described by 0 in Eq. (4); 0 (kp) = i/T+T^- At small T 
the function / is linear in T, and it saturates at T 1. 

The results are illustrated with the help of Eig. 3, 
where regions I, H and HI with a different dephasing 
rate behavior are indicated in the (T — p) plane. The 
regions are divided in accord with the importance of the 
renormalized spectrum of the FP and the relative contri¬ 
butions of the self-energy and vertex diagrams. In region 
I, which is on the left of the black line (i.e., at small 
densities), the characteristic momenta of p and q in Eq. 
(7) do not exceed qc- Therefore, the renormalization of 
the FP spectrum is important, and ojq ^ 52 - 77/2 
be used [32]. In region H, since the characteristic mo¬ 
menta of the FPs are larger than qc, it suffices to use 
the quadratic spectrum for FPs. In region HI, which is 
in the bottom part below the blue line, the dephasing 
time is long and only the self-energy diagram is impor¬ 
tant. Hence, the factor reduces to 1, and dephasing 
rate coincides with the out-scattering rate, obtained 
from the golden rule [5]. (In this calculation, qc just pro¬ 
vides an infrared cut-off.) Above the blue line, in regions 
I and H, both the self-energy and vertical diagrams are 
relevant, and the factor (t) is important; see also [33]. 
Due to the two-phonon structure of the correlation func¬ 
tion of the FP pairs participating in the inelastic process, 
the influence of this factor on the dephasing rate is rather 
non-trivial, so that one cannot expand (t). 

In Fig. 3, the blue and black lines have been found by 
matching the asymptotic behavior [28] of the dephasing 
rates deep in regions I, H and HI. We introduce (po, Tg), 
the values of the crossing point of the blue and black 
lines as characteristic scales: po and Tg ^ ^Po- 

Here, we have introduced 7 = ~ 0.02, which is a 

IP 
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parameter describing the adiabaticity of the el-FP inter¬ 
action. Under a given choice of parameters, it can be 
found numerically that ~ 0.02eV and Tq ~ 0.6K. The 
dephasing rate in different regions can be expressed as 

f 0.48 ifi/ no)^ I 

(T) ='yTx I O.lSyW^ II (12) 

[ 0.24^ logr' III. 

These expressions are obtained using asymptotic behav¬ 
ior of the function / in Eq. (11) and, therefore, are only 
applicable far away from the borderlines. At low enough 
temperatures, (t) = 1 and the function / (T, ^) is in¬ 
dependent of T. Hence, ~ in region III, which 
is a GR result. At high temperatures the phonons con¬ 
tributing to the electronic dephasing become quasi-static 
and, consequently, the dephasing rate is smaller than the 
out-scattering rate T~^^. Unlike region III, in regions I 
and II the dephasing rate is determined by a non-GR ex¬ 
pression, and is proportional to temperature, irrespective 
of rj. The existence of the linear in T regime is the main 
result of our paper. 

By comparing the rates in regions II and III, one may 
conclude that there should be a maximum in the dephas¬ 
ing rate as a function of fj,. Indeed, as it is illustrated 
by Fig. 4 such a maximum exists. The line indicating 
the maximum essentially overlaps with the borderline be¬ 
tween the regions I, II and the region III, which is illus¬ 
trated by the blue line in Fig. 3. 



0.05 0.15 0.25 0.35 

Figure 4. Dephasing rate as a function of the chemical poten¬ 
tial at different temperatures. From top to buttom: T=15K, 
lOK, 5K. 

Dephasing due to vector potential fluctuations. Un¬ 
like the scalar potential, the dephasing rates induced by 
vector potential are different for different channels ow¬ 
ing to the factor (f^,t) = 1 + where Si = —1 

for the intervalley Gooperons {I = 0, z) and s; = 1 for 
the intravalley Gooperons (/ = x,y). For the intervalley 
channels, the only relevant for the magnetoresistance at 
weak fields, the dephasing rate produced by the vector 
potential coupling is quite similar to its scalar counter¬ 
part, Eq. (12), with obvious modifications due to the 
change in the coupling constant and absence of screening 
for the vector potential [28]. 


Discussion. We have analysed the dephasing rate in¬ 
duced by FPs in graphene, and evaluated it for a wide 
range of n and T (see Fig. 3.) We determined several 
asymptotic regions with temperature dependence evolv¬ 
ing from ~ to ^ T when temperature in¬ 
creases. (See Fig. 4 in [28] for an illustration of the tem¬ 
perature dependence of the dephasing rate.) The tran¬ 
sition to linear behavior in T is related to the fact that 
at high temperatures phonons become slow on the time- 
scale of r^. 

The measured dephasing rate in graphene is usually 
compared to the contribution induced by the electron- 
electron interaction, , which is linear in T for T < 
1/Ttr [11]. However, the observed rate [12-14], when it 
is linear in T, always exceeds the theoretical estimation. 
In view of the linear dependence on T of the FP’s con¬ 
tribution to dephasing, it is reasonable to compare its 
value with rfl}. In principle, it is a competition between 
two mechanisms, each determined by a small parameter: 
the adiabatic parameter 7 and sheet resistance pn mea¬ 
sured in units of the quantum resistance. We compare 
the dephasing rates at density n = when the 

sheet resistance « Q.bkQ,. Under these conditions, both 
parameters 7 and pn are of the same value. Combin¬ 
ing the contributions arising from the scalar and vector 
potentials, we obtain Tflp/rflfl pc 0 . 2 . 

The in-plane phonons generate a dephasing that at 
T < is negligible compared with rflp, while at T > 
T^q the rate ~ T is comparable with rflp. (Note 
that for in-plane phonons, a region of non-GR dephasing 
rate, analogous to region H, develops at temperatures 
> p that is too high to be relevant.) It is important 
that each of the three rates Tpp, and has a 
distinct dependence on the chemical potential. While 
Tflfl decreases with density, rflp cc and oc p. 
This opens a way to identify each of these mechanisms 
by studying the magnetoresistance as a function of the 
chemical potential. 

In our consideration, we had in mind suspended 
graphene. However, our result may also be relevant for 
supported samples so long as they are coupled to the 
substrate by weak Van der Waals forces [34]. One may 
expect that such a weak coupling does not provide an 
essential change in the phonon spectrum. Indeed, it is 
known that the phonon spectrum in graphene [35] and 
graphite [36] are practically identical for the correspond¬ 
ing branches. FPs in supported samples have been dis¬ 
cussed recently in connection with the heat transport 
measurements in Refs. [ 6 , 7]. Until now flexural phonons 
have been a delicate object to detect in electronic trans¬ 
port. We propose here to observe them through weak- 
localization measurements. 

Acknowledgements. The authors gratefully acknowl¬ 
edge A. Dmitriev, I. Gornyi, V. Kachorovskii, D. Khmel- 
nitskii and A. Mirlin for the useful discussions and valu¬ 
able criticism. The authors thank the members of the In- 










5 


stitut fiir Theorie der Kondensierten Materie at KIT for 
their kind hospitality. A.F. is supported by the Alexan¬ 
der von Humboldt Foundation. The work is supported by 
the Paul and Tina Gardner fund for Weizmann-TAMU 
collaboration, and National Science Foundation grant 
NSF-DMR-100675. 

* K.Tikhonov and W.Zhao contributed equally to this 
work. 


* tikhonov@physics.tamu.edu 
^ wei.zhao@tamu.edu 

[1] S. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, 
Reviews of Modern Physics 83, 407 (2011). 

[2] A. H. Castro Neto, N. M. R. Peres, K. S. Novoselov, and 
A. K. Geim, Reviews of Modern Physics 81, 109 (2009). 

[3] E. Mariani and F. von Oppen, Phys. Rev. Lett. 100 , 
076801 (2008). 

[4] H. Ochoa, E. V. Castro, M. I. Katsnelson, and F. Guinea, 
Phys. Rev. B 83, 235416 (2011). 

[5] I. V. Gornyi, V. Y. Kachorovskii, and A. D. Mirlin, Phys. 
Rev. B 86, 165413 (2012). 

[6] A. A. Balandin, S. Ghosh, W. Bao, I. Galizo, D. Tewelde- 
brhan, F. Miao, and C. N. Lau, Nano Letters 8, 902 
(2008). 

[7] J. H. Seol et ai, Science 328, 213 (2010). 

[8] K. 1. Bolotin, K. J. Sikes, J. Hone, H. L. Stormer, and 
P. Kim, Phys. Rev. Lett. 101 , 096802 (2008). 

[9] E. V. Castro, H. Ochoa, M. 1. Katsnelson, R. V. Gor¬ 
bachev, D. G. Elias, K. S. Novoselov, A. K. Geim, and 
F. Guinea, Phys. Rev. Lett. 105 , 266601 (2010). 

[10] E. H. Hwang and S. Das Sarma, Phys. Rev. B 77, 115449 
(2008), note that the phonon’s contribution to resistivity 
is about 100 Ohm on the background of few kOhms. 

[11] B. L. Altshuler, A. G. Aronov, and D. E. Khmelnitsky, 
Journal of Physics C-Solid State Physics 15 , 7367 (1982). 

[12] X. Wu, X. Li, Z. Song, C. Berger, and W. A. de Heer, 
Phys. Rev. Lett. 98, 136801 (2007). 

[13] F. V. Tikhonenko, A. A. Kozikov, A. K. Savchenko, and 
R. V. Gorbachev, Phys. Rev. Lett. 103, 226801 (2009). 

[14] M. B. Lundeberg and J. A. Folk, Phys. Rev. Lett. 105 , 
146804 (2010). 

[15] J. Jobst, D. Waldmann, 1. V. Gornyi, A. D. Mirlin, and 
H. B. Weber, Phys. Rev. Lett. 108, 106601 (2012). 

[16] D. K. Efetov and P. Kim, Phys. Rev. Lett. 105, 256805 


( 2010 ). 

[17] J. C. W. Song, M. Y. Reizer, and L. S. Levitov, Phys. 
Rev. Lett. 109, 106602 (2012). 

[18] J. von Delft, F. Marquardt, R. A. Smith, and V. Ambe- 
gaokar, Phys. Rev. B 76, 195332 (2007). 

[19] P. M. Chaikin and T. C. Lubensky, Principles of con¬ 
densed matter physics, Vol. 1 (Cambridge Univ Press, 
2000 ). 

[20] H. Suzuura and T. Ando, Phys. Rev. B 65, 235412 

( 2002 ). 

[21] J. L. Manes, Phys. Rev. B 76, 045430 (2007). 

[22] D. R. Nelson and L. Peliti, Journal De Physique 48, 1085 
(1987). 

[23] P. Le Doussal and L. Radzihovsky, Phys. Rev. Lett. 69, 
1209 (1992). 

[24] K. V. Zakharchenko, R. Roldan, A. Fasolino, and M. 1. 
Katsnelson, Phys. Rev. B 82, 125435 (2010). 

[25] B. L. Altshuler, D. Khmelnitzkii, A. 1. Larkin, and P. A. 
Lee, Phys. Rev. B 22, 5142 (1980). 

[26] B. L. Altshuler, A. G. Aronov, and D. E. Khmelnitsky, 
Solid State Communications 39, 619 (1981). 

[27] E. McCann, K. Kechedzhi, V. I. Fal’ko, H. Suzuura, 
T. Ando, and B. L. Altshuler, Phys. Rev. Lett. 97, 
146805 (2006). 

[28] See Supplemental Material for the diagrammatic calcu¬ 
lation, evaluation and aymptotic properties of the decay 
function, which includes Refs. [37, 38]. 

[29] W. Eiler, Journal of Low Temperature Physics 56, 481 
(1984). 

[30] F. Marquardt, J. von Delft, R. A. Smith, and V. Ambe- 
gaokar, Phys. Rev. B 76, 195331 (2007). 

[31] V. N. Kotov, B. Uchoa, V. M. Pereira, F. Guinea, and 
A. H. Castro Neto, Reviews of Modern Physics 84, 1067 
( 2012 ). 

[32] We, however, neglect the h’^-vertex corrections. 

[33] G. Montambaux and E. Akkermans, Phys. Rev. Lett. 95, 
016403 (2005). 

[34] A. K. Geim and I. V. Grigorieva, Nature 499, 419 (2013). 

[35] D. L. Nika, E. P. Pokatilov, A. S. Askerov, and A. A. 
Balandin, Phys. Rev. B 79, 155413 (2009). 

[36] R. Al-Jishi and G. Dresselhaus, Phys. Rev. B 26, 4514 
(1982). 

[37] A. G. Aronov, A. D. Mirlin, and P. Wolfle, Physical 
Review B 49, 16609 (1994). 

[38] A. G. Aronov and P. Wolfle, Phys. Rev. B 50, 16574 
(1994). 



arXiv:1505.04834vl [cond-mat.mes-hall] 18 May 2015 


Supplemental Material for 

’’Dephasing time in graphene due to interaction with flexural phonons” 

Wei L.Z. Zhao/’^ Konstantin S. Tikhonov,^’and Alexander M. Finkel’stein^’® 

^Department of Physics & Astronomy, Texas A&M University, College Station, TX 77843-4242, USA 
^Department of Condensed Matter Physics, The Weizmann Institute of Science, 76100 Rehovot, Israel 
D. Landau Institute for Theoretical Physics, 117940 Moscow, Russia 
"^Moscow Institute of Physics and Technology, 141700 Moscow, Russia 
^Institut fiir Nanotechnolsoqie, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany 

(Dated: May 20, 2015) 


DIAGRAMMATIC CALCULATION OF THE DECAY FUNCTION 


The decay function F* defined in Eq. 
diagram: 


(6m) (’m’ refers to the main text) can be illustrated with the following 



FIG. 1: Diagrammatic representation of the decay function F (t). Here + EJ''^’’*. 


To the lowest order in the interaction propagators, the expression for F* can be obtained as [1] 




C{ 

cuti-t2y 


( 1 ) 


where C[ is the first order correction to the Cooperon in terms of the interactions propagators. We perform the 
diagrammatic calculation with the use of the matrix Greens function (cjP) = I have 

introduced iso-pseudospin basis [2], spanned by two sets of mutually commuting matrices, isospin-S and pseudospin- 
A: 


^0; Fix — G o'x^ Fy — G (Jy, — IIq G CT^, 

~ G (Xz, Ay = Ily G (Jz: A^ = Tlz G ctq, (2) 

where Eo,Ao are unit matrices and nx,y,z and (Jx,y,z are Pauli matrices acting on valley and sublattice spaces re¬ 
spectively. For each channel, the contribution of the el-FP interaction can be separated into the scalar and vector 
potential ones. Thus, one can write F* = F^ -|- F^. We concentrate first on dephasing caused by the scalar potential 
fluctuation (p (Q, Ut). The calculation for the vector potential contribution goes along the similar lines, and results are 
presented in the end of this Section. Due to the interaction propagators, one can write [1] 

c[ {h,h) = I {dQ) (dn) (dq) (dCj) [C' (q, Cj) {Co, Ul) C' {q, Co) 

+ {q,Co-n) E^7‘ ((h, fl) cl, {q,Co + n)], 


(3) 




2 


where ti 2 = — ^2 and ri 2 =^ 1 +^ 2 - Here, the effect of the interaction is separated into the self-energy and vertex 

(vertical) contributions, denoted as 2 . The self-energy and vertex contribution of the channel I are 

defined in the iso-pseudospin basis via y;bseH/t'ert = {Y]yh.yKi)^p {'^y^i^y)s'y (summation for Greek letters 

is implied; indices (^{A) are omitted). Here, the self-energy contribution, see Fig. 2(a), can be written as 


(a) 


(b) 


(c) 




=-i/2Di^{n,Q) I =l/7rmT 


FIG. 2: Diagrammatic representation of (a) the self-energy contribution and (b) the vertex contribution, (c) Diagrammatic 
dictionary for the various objects involved: ~ fJ-A iljl2 and = kp ± q/2\ uq is the density of states at Fermi energy, and 

gi (Q) is the coupling constant with screening included. 


{ [G« (p+, i?+) {p+ -Q.E+- n) (p+, £+)] {p-,E-)^, 

+G« {p+,E+)^^ [G^ (P-^E-) G^ {p- -Q,E-- ft) G^ {p-,E-)]^^}, (4) 

where E^ = M ± w/2, p^ = kp ± q/2. Note that arguments q and d) are related to the propagation of the Cooperon, 
while Q and H describe the interactions with FPs causing the dephasing. The vertex contribution, see Fig. 2(b), is 

(9: b!) = J {Sp) (Q, H)) X 

[G^ (p+,e+ - n) G^ + [G^ (p",e" + H) G^ (p" - (5) 

where = E^ ± H/2. Due to the softness of FPs, the integrals in Eq. (3) converge at small frequencies fl <^T. Note 
that the assumption that H ^ T depends crucially on the fact that T which has been checked a posteriori. 

Therefore, when calculating the effect of interaction on the Cooperon propagators, one may take only the classical 
(Keldysh) component of the interaction, (Q, O) = —2Z()) (Q, Q), and send n (wfc) in El (k) to its classical limit, 

If the full quantum problem is considered, one needs to use order to incorporate Pauli principle due the 

presence of other electrons [3]. This leads to the same result for T 3> but ensures that quantum fluctuations do 
not lead to dephasing at zero temperature. 

In Eq. (3), the Cooperon variables {q,Cj) are small comparing to the electronic scales determining the el-FP 
interaction process, i.e., q kp and uj ^ p. Therefore, one can drop out (g, w) dependences from the interaction 
propagators E. Eurthermore, as far as T ^ TbGi the typical momentum transfer Q ^ 2kp 1/vpT. Therefore, 
dephasing can also be calculated assuming that electron’s motion during the interaction event is ballistic. After a 
simple calculation, one concludes for scalar potential 


•spl,self _ spljVert 


(Q,ft) 


( 6 ) 
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where {Q, O) = ^ (Q, fi) {Q). Here </> (Q, fi) is defined in Eq. (7in), and (Q) is defined in Eq. (9m). Plugging 
these results back to Eq. (1) leads to Eq. ( 8 m). 

In the case of vector potential coupling, one calculates a diagram similar to that in Fig. 2, and gets 

where <5^ (Q) = 6^^ — jQ'^ and {Q, O) = VpA^^ (Q, H) Ba (Q) ■ Here, the vector potential correlation function 

Aj (Q, {g2/vFf J (d^p) (d^q) 

xHip)Hiq)Sp,^{n,Q), (8) 

where Sp^q (H, Q) = J2± d (fl ± Wp ± Wq) d (Q — p — q) and the summation includes four different processes of 
emission/absorption of two FPs by an electron. Also, the factor describing the ballistic motion of the electrons 

Ba (Q) = ^ (l - iQ/‘^kFf) d {2kF - Q) ■ (9) 

The corresponding decay function can be rendered as (compare with Eq. ( 8 m)) 

Fa (^) = ^ / (Q’ iQ) Ba (Q) Cf (H, t) , (10) 


where 




1 + S; 


sin nt 


( 11 ) 


The intervalley Cooperons are coupled to the vector potential field of the opposite signs. Thus, s; = ±1 for intra- 
/intervalley Cooperons. Further on, one can resolve the transverse delta function dp and get 


F’a {t)=t J sin^ Q (dQ) (dO) A (Q, fl) Ba (Q) Cf (fl, t ), 


( 12 ) 


where 

A{Q,n) = ^J (d^p) {d\)p\^H{p)H{cDSp,q{n,Q)- 


(13) 


EVALUATION OF THE DECAY FUNCTION 
Scalar potential coupling 

Let us now evaluate the integral in Eq. (8m) explicitly. It is convenient to use the time representation for the 
energy delta-function: 


(27r) d (O ± Wp ± Wq) 


dr cos (ujpT) cos (w^r) exp [—iflr] . 


After this, one can integrate Eq. (8m) in frequency fl, using 


I 


(dfl) 



sin nt \ 

nt ) 



{r/t), 


where 

(a) ^d (a) ±20(1-1^1). 


(14) 


(15) 


( 16 ) 
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The next step is to make the integral dimensionless by introducing t = st and dimensionless 2D vectors 
X, y, z = p/feiT-, q/ZcF, Q/^f- Using the expression for H (q) in the classical limit, H (q) = {ojq is defined in 
Eq. (4m)), we obtain the decay function as given in Eq. (11m) 


(t) = cltTf^ (^akpt, Z ^^'^qc/kp 


r /O- A 2 f tn \ \ \ X t ^ ^ ^ COS {T sx^Q^ {x)) COS {T sy^Q^ (y)) 

U 0 = 47 r J (2 f) (dz) (dx) (dy) d (z - x - y) S' (z) J (s) ds -^- 0 ^ (J) 0^ (y) -^ 


with T = ak%t^ ^ ^ ^^^Qclkpj (a:) = y 1 H- (x/^) and 


S 0 (z) = 


+ zj (geN) 


'l-izl2fei2-z), 


which is the product of the screening and chiral factors. Here = 4 is the spin-valley degeneracy in graphene, and 
ge describes the renormalized Coulomb interaction. 

To proceed, it is convenient to make some transformations in Eq. (18). First, we integrate out s exactly, using the 
relation 


J (s) ds cos (so) cos (sb) = S± (a, 6 ), 


E:± [a, 5] = 1 ± 


i sin a cos b — b sin b cos a 

a? — 52 


This gives 


t tn- ^4 //J ^ w , , , (x-y)^^- [T05 (x)a:2s,r0^ (?/)?/ 2 s 

{T,0 = ( 2 ^) J (dz) (dx) (dy) d (z - X - y) S,^ (z) - q 2 q 2 - 


The final step is to resolve the delta-function in the above expression for y which yields 

y = \/z'^ + x"^ — 2zx cos tp, 

where ip is the angle between z and x. As a result, one finally obtains 

tn- ^ c ( \ ^ dp;{zcos'ip-xf'E._[rQi{x)x‘^,rQ^{y)y‘^] 

U(T.t) = -efMefM- 


Vector potential coupling 


We define 


F\ {t) = c\tTf\ (aklt,Z-^/'^qJkp) -, 


where ca is the dimensionless el-FP coupling constant for the vector potential. 


Intervalley channels {I = 0, z) 


Since si = — 1, the decay function for the intervalley Cooperons is similar with that for the scalar potential. Without 
providing further details, we conclude that 


fA{r,0= / dzSAiz) dx 


/■“, dp, 1 S_ [r0c(a:)x^r0;(d)d^] 

0 ^ Jo 2Fxy2 02(a;)0|(y) 
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where y has been defined by Eq. (23), and 

SA{z) = l{l-{z/2fy'^%{2-z). (27) 

Note that Sa (z) includes only the chiral factor since screening does not affect the vector potential coupling constant. 


Intravalley channels {I = x, y) 


For the intravalley channels, s/ = 1. Correspondingly has to be changed to 


c 2 -k 


fAir,o= / dzSAiz) 


dx 


dip 1 :z+[T<d^{x)x‘^,T&i{y)y^ 


27r xy^ 


02 (x) 02 (y) 


(28) 


ASYMPTOTIC PROPERTIES OF THE DECAY FUNCTION AND THE PHASE DIAGRAM 

The analytical expressions obtained in the previous sections allow us to determine t^/a arbitrary tempera¬ 
ture T and chemical potential y. (In the end, the dephasing times t^/a are defined as solution of the equations 
d^ip/A {ftp/ a) = !■) We have developed numerical procedure, which exploits these equations to calculate correspond¬ 
ing (T, /i) dependencies. Before presenting general results, let us concentrate on the properties of the functions 
f 4 >/A (T,^) in the analytically accessible regimes. 


Scalar potential coupling 


Large ^ 


For ^ 1, one may identify three asymptotic regions depending on T. When T<C C the integral in (T, $) is 

dominated by the quadratic spectrum where 0^ (x) = 1 (refer to Eq. (24)), while at <C T <C the integral 

is dominated by the part of spectrum where 0^ (x) = {x/^)~^^^. Thus, for not too large T {T one may 

assume that the spectrum is homogeneous, and put 0j {x) = ix/ppf’. Then the function reads: 


U{r,o = 


4b 


27r 


/ dzS^ (z) / 

/O JO 


Jo 


\-\-2b yA-\-2b 


X 


T{x/Ox\r{yKfy" 


We obtain, 


6b 2+45 

/0(r,e) = c(&)^5Trr^, 


(29) 


(30) 


with c ( 6 ) = ^ 2 +b ’ where dzS^ (z) = 0.89 (assuming ge « 1), and 


Mb 


_ 5b+4 

U ^ + 2 ■ 


!_ {u, u)du = 2 *>+2 r 


55-|-4\ . (56-I-4) TT 

b + 2 ) {b + 2)2' 


(31) 


More specihcally, for 6 = —? 7 / 2 , Cj = c^—yjlP) = 0.03, when r] = 0.8. For 6 = 0, C// = c(0) = fiV^. The scaling 
expression in Eq. (30) can be achieved due to the fact that one can neglect all the z dependences except for (z) 
in Eq. (29). Then, y "Zi x, ip integral gives 27r, and the remaining integral over x leads to Eq. (30). The obtained 
asymptotes will describe non-golden rule behavior for the decay function, in the region I (6 = —? 7 / 2 ) and II (6 = 0). 
Note that for 6 = 0, /^ = c//T. 

Next, for T the golden rule regime holds and the vertex diagram is not important any more. In this case, 

one can put S_ -A I, and (T, C) becomes independent on T: 


U{r,o 




dip {z cos ip — x)^ I 
2tt xy^ 0 | (x) 0 | {y) 


( 32 ) 
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where c^y = t]U^W^ dzS^ (z) Here, with J„ (u) being the Bessel function of the first kind, = 

u-‘^^Jiiu)du = J2 (u) du)^ = + 


^2’)-ir(l+7;/2) V 


r{2-v/2) J 


For r] = 0.8, this yields = 1.63 and = 2.42. Eq. (32) describes the asymptotic behav¬ 
ior of the function (T,^) in region IV. 

As a result, we see that the function (T, for large ^ evolves with growing T as follows (regions of applicability 
can be easily read off from these equations): 


ciiT ^ ^ (33) 

As indicated in the main text, the matching of different asymptotes gives us a separation between different regions of 
the dephasing rate on the (T — /i) plane: The black line in the Fig. 3m corresponds to matching of the first pair of 
asymptotes. The location of the region IV is related to matching of the second pair. 


Small 5 

For small ^ 1, there are two asymptotic regions depending on the value of T. For T ^ 1, the integral in Eq. 

(29) is always dominated by the quadratic spectrum (& = 0), and function (T, C) acquires the asymptote linear in 
T (compare with Eq. (30) at 6 = 0): 

UiT,0 = cnT. (34) 

At large T ^ 1, the golden rule is applicable. In this region, 

f<j> {T, 0 = cm log 1/C, (35) 

where Cm = (z) jz^dz. Note, that Cj, c//, Cm and Cjy implicitly depend on via Stf, (z), see Eq. (19). Thus, 

in the case of small C, (T, C) evolves with growing T as follows 

C 77 T Cm log 1/C- (36) 

The blue line in the Fig. 3m corresponds to matching of the above asymptotes. 

Equation (24) allows for straightforward numerical evaluation. Here as an illustration, we present the numerical 
calculation of the function (T). From the plot we see an excellent agreement with theoretical calculation for the 
region T <C 1 and T 1, even for an intermediate value of C = 0.52. 


h 



FIG. 3: Plot of the function (T). The parameters are taken as follows: rj — 0.8, C = 0.52, and Z = 3. The blue dots are 
numerical calculation and the red lines are the asymptotes for different region mentioned above. 


Phase diagram of the dephasmg rate 

In Fig. 3m, the blue and black lines have been found by matching the asymptotic behavior of the dephasing rates 
deep in regions I, H and HI. The blue line is found as a consequence of matching Eq. (34) with Eq. (35). It is defined 
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by the criterion uj^pT^ ~ 3, where is calculated numerically from Eq. (11m) with the phonon frequency tOkp taken 
from Eq. (4m). The black line is found by matching the first pair in Eq. (33), and it can be obtained according to 
the equation ~ 7. As mentioned in the main text, the obtained rate has a transition from ^ to ~ T 

when temperature is increased. This transition is illustrated on Fig. 4. 


T0 ' (ps ') 



FIG. 4: Depahsing rate as a function of temperature at two different chemical potentials. The transition from a quadratic to 
linear dependence takes place at ~4K for the red solid line. For lower chemical potential (represented by the blue dashed line), 
the transition happens at a lower temperature as shown in the inset. 

Note that besides the regions discussed in the main text, on the {T — /i) plane there exists another asymptotic 
region (IV), not mentioned in Eqs. (12m). It lies at very small densities below the blue line and above the line ^ = I. 
In this region, the dephasing rate is still described by the GR, but contrary to region III, renormalization of the FP 
spectrum at low momenta is essential. In this situation, the dephasing rate equals ~ qT ^ (T/Tq)^ ^. 

However, the relevant densities are so small, that this region does not fit into the scale of Fig. 3m. 

It is instructive to compare the situation with flexural phonons to that for ordinary phonons in semiconductors, 
where the dephasing process can be viewed as coming from the energy diffusion via low-energy (i.e., quasi-elastic) 
collisions [4]. Under these conditions, the dephasing rate \/ (j<i>lTout) where 5e denotes the characteristic 

energy transfer during a single scattering event. Thus, the accumulation of phase in the course of the energy diffusion 
yields ~ /Tout)^^^■ The specific point of our problem, compared to dephasing via usual phonons, is that 

an electron is coupled to two flexural phonons, and dephasing cannot be described by the energy diffusion process. 
Because of this two-phonon interaction, the support of the correlation function of the fluctuations is not characterized 
by any typical frequency and frequency transfer occurs in such a way that typical energy transfer is of the order of 
r^^(T). In particular, this allows for an energy transfer exceeding TbGi although the momentum transfer is limited 
by 2kp. As a result, the decay function F^{t) at short times (i.e., in the non-GR regime) is proportional to rather 
than as for the case of the energy diffusion. Note that the energy diffusion corresponds to the expansion of the 
factor C'^(<) ^ (Hf)^ which leads to the t^-dependence of F^{t). As we have already mentioned, in the case of FPs 
one cannot expand C^{t). 


Vector potential coupling 

Intervalley channel 

The similarity (compare Eq. (24) with Eq. (26)) with the scalar potential case makes the calculation of asymptotes 
for the function (T, C) straightforward. The function (T, ^) evolves with the growth of T in the following fashion: 
for ^ » 1, 


diiT ^ ^ (37) 

while for f ^ 1, 

diiT—t diii^~^. (38) 

Note that due to the absence of screening, dependence on ^ is stronger than in the case of the scalar potential, compare 
Eq. (36) with Eq. (38). 








Let us present some technical details. Here dJj = d{—ri/2 ), where d{b) = and is defined in Eq. (31); 

specifically, for 6 = 0, du = d{0) = 1/64. In the GR region, the function (T, C) is independent on T. Its asymptotic 
behavior is given as follows: 


fAiO 


dzSA (z) 



dip 1 1 

27r xy'^ 0| (x) 0| (y) 


f diiiP, 1 for ^ < 1 
\ d]yp-^^ for C » 1 ■ 


Here, for r; > 1/2 (which is important for the convergence of the following integral). 


, 1 , /■“ , dip 1 1 

2 X ^ X ^ X 27r xy 2 02 ^^ [x) 0 |^^ (y) ’ 


and 


d]y = 2yUW^GA, 


(39) 


(40) 


(41) 


where is the same as in the case of the scalar potential, ^ Jo {u) du)^ 

dzSA (z) d — yields = 1.68 and Ga = 1.73. 


/ 2^-H(r,/2) 

V r(l-7)/2) 


2 

and Ga = 


Intravalley channel 

For the intravalley Cooperon, the function {T,0 involves which does not vanish in the limit of T —> 0. 

Therefore, it evolves with the growth of T (i.e., from T <C 1 to T 1) as 


2/a (?) ^ /a (?), 


(42) 


where pA (?) is defined in Eq. (39). It is noteworthy that the static deformations produce dephasing here. This 
happens because the electrons on the interfering trajectories are coupled to the vector potential of the same sign. The 
situation is similar to that for a particle in the random magnetic field [5, 6 ]. Unfortunately, the observation of this 
effect is obscured by the gaps inevitably produced in these channels by disorder scattering. The rate obtained from 
(^) = 1 changes from 2 t“/^ to t~\ with increasing temperature, where 


' *,A 


= 7T 


for?»l 

1.9y/T/^ for ? < 1 


(43) 


Here 7 ~ 0.02 is the adiabatic parameter defined in the main text, /Tq ~ jAc and T/ ~ 7 ^Ac where Ac describes the 
energy scale of the anharmonicity (see the main text). Note that for ? <C 1, the characteristic momentum transfer 
is ~ Qc instead of kp. As a result, the rate is independent of the chemical potential. Note also that for 

the intravalley channels the diffusive limit {qj 1) for the calculation of dephasing due to the el-FP interaction can 
be achieved at higher temperature as compared to the intervalley channels. For example, for dimensionless sheet 
conductance ~ 10 and chemical potential y ~ O.leV, the diffusive limit occurs already below T ~ 0.5K. 
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